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Abstract- Representative atomic and molecular systems, including various inorganic and organic molecules with covalent and ionic bonds, have been studied by 

using density functional theory. The calculations were done with the commonly used exchange-correlation functional B3LYP followed by a comprehensive analysis of 
the calculated highest-occupied and lowest-unoccupied Kohn-Sham orbital (HOMO and LUMO) energies. DFT calculations of vibration spectra of many organic 
systems [1, 2], Density Functional Theory (DFT) calculations have been performed to support our wave number assignments. The theoretically predicted IR is well in 
agreement with that of experimental spectral data. 

♦ 

1. Introduction 

The global properties may explain the reactivity and selectivity. One of these global quantities is hardness g. The relation of the hardness with the HOMO- 
LUMO gap is physically clear. It is obvious that the more stable structure has the largest HOMO-LUMO energy gap. Therefore, an electronic system with a larger 
HOMO-LUMO gap should be less reactive than one having a smaller gap. It is the base of the maximum hardness principle. Another global property is the electric 
dipole polarizability. 

Hardness measures the stability and polarizability measures the reactivity [8-9]. The minimum polarizability principle (MPP) [15] has been postulated which 
expected that 'the natural direction of evolution of any system is towards a state of minimum polarizability'. Molecular properties such as dipole polarizabilities have 
been studied in thiophene monomers and polythiophenes by many methods . 

2. Methods 

2.1 Computational method The ionization potential is calculated S — 1/2tj 

in present work as the energy difference between the energy of the =((d A 2 N)f l3Ej A 2 )_(" v" (r)) = (3N/&k)_C v " CO) ■» ( 2 - 9 !) 



molecule derived from electron-transfer (radical caption) and the 
respective neutral molecule [2, 3, 4, 5, and 6]: 
*P = ^^-^..(2.1) 

The EA was computed as the energy difference between the 
neutral molecule and the anion molecule: 
ujH — £ n ^cnLon t i,i J 

The HOMO and LUMO energy was also used to estimate the IP 
and EA in the framework of Koopmans' theorem [6, 7, 8, 9, and 
10]: 

IP = -Z-H9M9 — ( 2l3 !l 

" = -Zlumo 

2.2 Chemical potential , chemical hardness and softness 

Within the framework of the density functional theory (DFT) , 
one of the global quantities is chemical potential (k), which 
measures the escaping tendency of an electronic cloud .It is a 
constant , through all space , for the ground state of an atom, 
molecule or solid , and equals the slope of the energy versus N 
(number of electrons) curve at external potential v(r) [9,6,10,11]: 
jc = [&EfdN\_C f pr) 3 — CZ- 5] 

Finite difference approximation to Chemical Potential gives, 
pc ns -x = -{{IP + EAJ)/2... (2.6) 

Where % is electro negativity, the electro negativity is a measure of 
the tendency to attract electrons by an atom in a chemical bond 
and defined as the negative of the chemical potential in DFT [7,18] 
The chemical hardness is a measure of the resistance to charge 
transfer. The theoretical definition of chemical hardness has been 
provided by the density functional theory as the second derivative 

Finite difference approximation to Chemical hardness gives, 
17 = QP-EA)/2«. (2.8) 5 

The chemical hardness of a given system is , as can be seen by 
equation (3.13) , a resistance of the chemical potential to change in 
its number of electrons . For insulator and semiconductor , 
hardness is half of the energy gap. 

The global chemical softness, S, is a property of molecules that 
measures the extent of chemical reactivity.lt is the inverse of the 
chemical hardness rj [12,13,14]: 



2.3 Electrophilic Index 

Electrophilicity is a measure of energy lowering due to maximal 

electron flow between donor and acceptor. Electrophilicity (go) is 

defined as [13,15,16, 17]: 

qj = K A 2/(2r? ) ... (2.10) 

3. Calculation Method 

Geometries of all of the molecules considered in this study were 
fully optimized by using gradient corrected DFT with Becke's 
three-parameter hybrid exchange functional and the Lee -Yang- 
Parr correlation functional (B3LYP)t 19 ]we also performed single 
point energy calculations on the neutral systems using basis sets , 
Semi empirical methods Pm6, Hartree Fock with basis set 6-31 1G* 
and density functional theory the B3LYP with a larger basis set , 6- 
311G* 6-311G**, aug-cc-pVDZ , aug-cc-pVTZ,aug-cc-pVQZ , SDD 

4. Results and Discussion 

Basis Set Dependence of the DFT Results . Before exploring any 
correlation of the DFT HOMO /LUMO energies with molecular 
properties , we need to understand the basis set dependence of the 
HOMO/LUMO energies . Previous DFT calculations 15 tested a 
series of correlation-consistent basis sets , starting from the 6- 
31 1G*, 6-311G**,aug-cc-pVDZ,aug-cc-pVTZ, aug - cc - pVQZ , 
SDD to obtain HOMO and LUMO energies for developing 
correlations with molecular properties . 

We note that in the previous study, 4 quite small molecules were 
used and there was no issue in terms of the computational 
resources needed, whereas we are interested in developing 
correlations for much larger molecular systems where a smaller 
basis set is needed in order to perform efficient computations. 

The calculated results are collected in Tables together with 
available experimental data. Note that the calculated vertical IP 
and EA values listed in the tables 1 and 2 were determined from 
the total energy calculations on the neutral systems. 
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Table (1): Experimental and Calculated Ionization Potential (IP) Table (6): Calculated softness for molecules 

for molecules 
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Table (2): Calculated energy lumo invertion for molecules 
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Table (3): Calculated energy gap for molecules 
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Table (4): Calculated chemical potential for molecules 
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Table (5): CalculatedChemical hardness for molecules 
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Table (7): ElectrophilicityCalculated for molecules 
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Conclusions 

Energy gap and hardness presented an assistance to predict 
reactive molecules. 

The ,Semi empirical methods Pm6, HartreeFockpresented good 
results in calculation oftotal energy 

he ionization potentialand electron affinity of a molecules 
calculated by using energy-vertical method werefound in a good 
agreement with experimental results and better than those 
calculated by using KT/DFT and KT/AM1, thus Koopman's 
theorem was not satisfied accurately. 
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